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Abstract 

A standard Lax representation for the polytropic gas dynamics is de- 
rived by exploiting various properties of the Lucas and Fibonacci polyno- 
mials. The two infinite sets of conserved charges are derived from this rep- 
resentation and shown to coincide with the ones derived from the known 
non-standard representation. The same Lax function is shown to also 
give the standard Lax description for the elastic medium equations. In 
addition, some results on possible dispersive extensions of such models are 
presented. 



*email: alexcOpas. rochester.edu 
temail: das@pas.roch6ster.edu 
■l-email: toppanacbpf .br 



1 



1 Introduction: 



Systems of hydrodynamic type [l]-[5] have been studied extensively over the last 
several years. Such systems manifest, among others, in string theories, mem- 
brane theories and topological field theories [6], The polytropic gas equations 
[7] belong to this general class of systems and are described by 

ut + uux + v'^~'^Vx = 0, vt + {uv)x = (1) 

where u, v denote the two dynamical variables. For different values of the expo- 
nent, 7, these equations describe different physical systems, which are disper- 
sionless and integrable [8]. 

Being dispersionless integrable systems, hydrodynamic systems can be de- 
scribed by a Lax equation on the classical phase space [8] . The Lax description 
that has been obtained so far for the polytropic gas equations [9]-[10], how- 
ever, is what is called a non-standard representation, which is not very useful in 
generalizing this system to other cases, such as the supersymmetric polytropic 
gas equations [11]. Nevertheless, the non-standard Lax description has been 
quite useful [9]-[10]. It has led naturally to the two infinite sets of conserved 
charges of the system (by the standard construction, although the two sets are 
obtained by expanding the residues at different points). The Lax description 
also immediately leads to the involution of the charges and clarifies why both 
the polytropic gas equations and the elastic medium equations share the same 
set of conserved charges (basically because they are both described by the same 
Lax function). One can also construct the Hamiltonian structures from this 
description, although a more convenient construction is through the Moyal-Lax 
representation of the system [12]. 

Many integrable systems can have both a standard as well as a non-standard 
Lax representation and, as we have already mentioned, a standard representa- 
tion is much more useful. Nonetheless, such a description is lacking so far and, 
in this paper, we construct such a representation for the polytropic gas equa- 
tions. Interestingly, such a description involves the use of Fibonacci and Lucas 
polynomials and their properties [13]. Although a standard Lax representation 
has unique residues, we show that such a Lax description, nevertheless, leads 
to the two infinite sets of conserved charges coming from two distinct families 
of fractional powers of the Lax function. The same Lax function also provides 
a standard Lax representation for the elastic medium equations [7, 9, 10]. 

The paper is organized as follows. In section 2, we review, very briefly, the 
definitions and some essential properties of the Fibonacci and Lucas polynomials 
[13]. In section 3, we derive some identities satisfied by two auxiliary functions 
dependent on these polynomials. Using these, we construct, in section 4, the 
standard Lax representation for the polytropic gas equations. In section 5, 
we construct the two infinite sets of conserved charges for this system, which 
coincide with the earlier known results [9]. The involution of these charges is 
automatic, since they come from a Lax description. In section 6, we attempt 
to construct some dispersive equations, whose dispersionless limit may lead to 
these equations. This question is extremely difficult and we present only some 
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partial results on this issue. We end with a brief conclusion in section 7. In 
the appendix, we show how the same Lax function that leads to a standard Lax 
description for the polytropic gas equations also describes the elastic medium 
equations. This also shows that the two systems share the two infinite sets of 
conserved charges. 



2 Definition and properties of the Lucas and Fi- 
bonacci polynomials: 

The Lucas polynomials are defined recursively as follows: 

ln+l{x) = X ln{x) + ln-l{x) (2) 

with lo{x) = 2 and li{x) = x. Their explicit form for n > 1 is: 

where \_x\ is the largest integer smaller than or equal to x and ^ ^ ^ is the 

binomial coefficient. 

The Fibonacci polynomials are defined by the same recursion relation: 

fn+l{x) =xfn{x) + fn-l{x) (4) 

but with fo{x) = and fi{x) = 1. Their explicit form for n > 1 is: 

k 



Ux) = J2 

k=0 ^ 



n-k-1 \ ^n-2fe-i 



A good presentation of these two families of polynomials and the many 
relations involving them can be found in [13]. However, for convenience, we list 
below some of the relations which we are going to use throughout the paper: 

Inix) = fn+l{x) + fn-l{x) (6) 
{X^ + 4) fn{x) = ln+l(x) + In-l(x) (7) 

ll{x)-ix^+i)mx)^A{-ir (8) 

l'n{x)=nfn{x) (9) 
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3 Various useful identities: 

In addition to the Lucas and Fibonacci polynomials, the following two families 
of polynomials are also useful in deriving some of our results: 



I — n ^ 



fe=0 

LfJ 



n — k 
k 



(10) 



m:^) = E("/)^'= (11) 



If a: ^ 0, one can easily see from (3) and (5) that: 

ln{x) = x^gn (12) 

for n > and: 

Ux) = x'^-'K-i (^^) (13) 

for n > 1 . Using this connection with the Lucas and Fibonacci polynomials, we 

can prove that 

Proposition 1 For any n>2, the polynomials Qn and hn satisfy the following 
relations: 

(i) hn{z) = hn-l{z) + Z hn-2{z) 

(ii) g„{z) = 5„-i(z) + z gn^2iz) 

(Hi) gn{z) = hn{z) + Zhn-2{z) = hn-l{z) + 2zhn-2{z) 

M 9niz) =nhn-2{z) 

(v) gn+i{z) - z g'„+-^{z) = gn{z) -^zg'^{z) 

(vi) ^,9niz)g'^^^{z) - l5;(^)5„+i(^) = (-^)"-^ 
Proof. To prove (i), we use (13) to rewrite (4) as 



for X ^ 0. If we now let z = x~'^ we find that (i) holds for z ^ 0. Since both 

sides of (i) are polynomial, it follows by continuity that it must also hold at 
z = 0. In a similar way (ii) follows from (2), (iii) follows from (6) and (iv) 
follows from (9). Let us note from (iii) and (iv) that 

2 

9n+l{z) - :^J-:pY^fln+l(-2) = hn{z) + 2 zhn-l{z) - 2 zhn-\{z) = hn{z) 
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as well as 



9niz) Zg'^(z) = hn{z) + Zhn-2{z) - zhn-2{z) = hn{z) 

n 

so that (v) holds. Finally, for (vi) we proceed by induction. For n = 2 the 
relation obviously holds. Now suppose that it holds for some n and notice that: 

= g-n+liz) hn{z) - hn-l{z) gn+2{z) 

= gn+i{z) {hn-i{z) + zhn-2{z)) - hn-i{z) {g„+i{z) + zgn{z)) 

= {-Z) {hn-l{z) gn{z) - gn + l{z) hn-2{z)) 

= \ 9n{z) - I- 9n+i{z) g'^{z) 

The relation therefore holds for n + 1, which completes the proof. ■ 
Let (f){z) be a function which is analytic in a disc D = {z : \z\ < r} and let 



n=0 



be its Taylor expansion around z = 0. We will denote by (2;) the polynomial 
consisting of terms in the Taylor expansion of (f> around 2; = up to z™. Let 
P{z) = ao + ai 2 + . . . + a„ be a polynomial. We will denote by P the polyno- 
mial P{z) = a„ + a„_i z + . . . + aoz"^ and we define P_^ (z) = z™' [-P^] ^ {z~^) 
for any positive integer m. With these notations in place, we can now state the 
following: 



Proposition 2 The following two relations hold for any n>2: 

{z)=gn+i{z) 

L 2 J 



(i) 
(ii) 



- 71 + 1 - 


n — 1 


- 71+1 - 

5"" 


- 71-1 - 


(-) = 

n—2 


- 71-1 - 

5™" 



,71-1. W =5n-l(2;) 

L 2 J 



Proof. For (i) notice that 5„" solves the differential equation 

J—g^{z)<l>'{z)-^g'„{z)<t>{z) = Q 

with initial condition 0(0) = 1. From (vi) of Proposition 1 we see that gn+i 
solves the differential equation 



n + 



^ -^gn{z)<l>'{z)-^g'^{z)<j>{z) = {-zy 



(14) 



5 



n+1 

with initial condition (f>{0) = 1. It follows from this that ip = g-n+i — gn" 
solves (14) with initial condition ^(0) = ip{0) = 0. If we now evaluate both 
sides in (14) at 2; = 0, with the identification (j) = tp, we get that V'(0) = 0- 
Differentiating (14) and evaluating at z = (again with the substitution ^ = ^) 
we get ^"(0) = 0. Repeating the process n — 2 times we get 

V)(o) = v^^Ho) = . . . = V^""'^(0) = 



It follows from this that 



yn+i 



Since Qn+i has degree [^^J and [^^J < n — 1 when n > 2, it follows that 



n-l 



(^) 



The proof for (ii) is similar. 



Proposition 3 For any integers n > 3 and k > let rn,k be the coefficient 
of in the Taylor expansion of ((7„(z^))'^+~ around z = and fn,k the 

coefficient of z''" in the expansion of {gn{z'^))''~ ^ ■ Then, 



rn,k = < 



n (mn+1), k = 2l 



I 0, 



k = 2l + l 



fn,k = < 



(-1)" 



I 0, 



n (mn-1), k = 2l 
m=;+i 

k = 2l + l 



Proof. We are only going to prove the formula for rn,k , since the proof for f„_fe 
is almost identical. 

For fc = 0, it is easy to see (using (iv) of Proposition 1) that 

d / — \ 1 / ——1 \ 

r„,o = \ 9n {z)j _^ = - n/i„_2(0) [g^ (0)j = 1 

For k = 1, r„,i is the coefficient of 2:"+^ in the Taylor expansion of g„ " (z^). 
If n is odd this is obviously 0, because the Taylor series contains only even 
powers. If n is even, namely n = 2m, then r2m,i is equal to the coefficient of 
^m+i ^jjg Taylor series of g2m^'" (-^)' which is according to Proposition 2, if 
n > 3. Hence, for any n > 3 , r„,i = 0. 
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Prom the residue theorem and the fact that gn{z'^) = z'^ln{z ^) it follows 
that: 



= jr^ira ^"(^") (-"«"(--^))'"^" (15) 

where the contour of integration F is contained in a neighborhood of ^; = 
where {gn{z^))''~^"^ is analytic. Prom (15) we obtain after integrating by parts: 



kn + 2 2iTi dz \^0fe"+2^ 

2kn+2 If 1 _i.Nfc-i+i 
= 2^i Jr ^(^^ 

where we have also used (6) and (4). We can rewrite this last result as: 
kn + 2 1 /" 1 J. /_ni 



2kn + 2 ' 27riir 
Moreover, it also follows from (15) that 

1 / 1 7 (.n, . -IxNfc-l+i 



= ^ / ;(^^«"(--^) (-"'n(.-^))* 
= 7(^3^ (/n+l(.-^) + /„-l(.-^)) {z-ln{z-'f-'^" dz 

''^hl ^04^^"+^(^") (-"^n(.-^))'=-^^" dz 



kn+2 

■Tr, 



2kn + 2 
and, therefore, we have 

^ V. = ^ / ^(^/n.i(^-^) (^"^4^-))'=-^^" (IT) 
Now if we multiply (16) by n — 1 and (17) by n + 1 and add them we get: 

X (0"Z„(-2-^))*'"^+^ dz 

= h / ^(-^i; (^-^(^")+^"-(^")) (^"^4^-))'=-^"^ 
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= _((;, _1)„ + 1)_L _L_ (,-2 + 4) 

_((fc_l),+ l)_L|^ _L_ (,-2 + 4) (,"i„(,-i))'=-2+* dz 



so that 

fcn^ + n — 1 1 /" 1 



(A;n+l)((A;-l)n+l) 2m z(^-^>'+^ 

= X {z''ln{z-^)f~^^" dz (18) 
Finally, adding (17) and (18) we get: 

4(-l)"r„,._, = ±-J^-^^^[Uz-^)-{z-'+4)f^{z-')] 

X (z"l„(z-i))''"^+^ d2 
kn {kn + 2) 



(fcn + 1) ((A; - l)n + 1) 



In other words: 



'^"•'^"^ > kn{kn+2) ^^^^ 

and the desired result follows by induction. ■ 

4 The Lax function and the Lax equation: 

In this section we show that, for n > 3, the Lax function 



provides a standard Lax description for the polytropic gas with parameter 7 = 
n+1, namely 

d Ln n 



-^{(4^^)^,^ (21) 



dt n+: 

where ()_|_ denotes the part of the polynomial that contains only non-negative 
powers of the variable, leads to the dynamical equations of motion 

Ut+uu^+v'^-'^v^ = (22) 
vt + {uv)^ = (23) 
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which are the polytropic gas equations. 

First, let us introduce the following simplifying notations: 

w = -V, 7r„ = p, Ln = Ln-u 

= e ^ - 1^ If = {ni «..)~^ {A,B] 

where {A,B} denotes the conventional Poisson bracket of A,B. With these 
notations, we can now rewrite (20) as 



Ln = - {^/wY In 



W TTj 



u=^Tr^gn{wTT„'^)+u (24) 



From this we see that 



1 + :^ 



1+^ 



n+1 



n+1 



{gn{w z^)+nu : 



J n+1 



TT, 



n+1 



{gn{w z^)f^ 

{9n{z)f^ 



n+1 



(^ = 7r-i) + (n + l)n7r-" 



(^ = «;7r„^) + (n + l)u7r„" 

L 2 J 



^n+l 



T^n" 9n+l (w7r„^) + (n+ 1) 



(ei)^ + („ + i) (i 



WKr, 



n+1 



(25) 



where we have used Proposition 2. Using the above relation, the Poisson bracket 
on the right hand side of the Lax equation (21) becomes: 



1 Ln 



d_ 

dp 



+ 



n+1 



UxP- 



dLn dL„ 



dp dx 



(26) 



Now, with (v) of Proposition 1 we obtain: 

■< i^n^A^T^n^) - («;7r-2)5;_^i(w;7r-2) 



n+1 



n 
n+1 
2n 
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n+1 ( 5L„ - \ 

= -Z TTn ^ \-nLn\ 

2n y diTn J 

n + 1 / aZ„ ~\ 

(27) 

Furthermore, using (vi) of the same proposition we notice that: 



Putting together (26), (27) and (28) we see that the Lax equation (21) is equiv- 
alent to 

Equating coefficients of powers of p in (29) we see that it holds if and only if 
the equations of motion (22) and (23) hold. 

This shows that the Lax function in (20) does provide a standard Lax de- 
scription for the polytropic gas dynamics. The same Lax function also provides 
a standard Lax description for the elastic medium equations, which we show in 
the appendix. For 7 = 3 (n = 2), the polytropic gas equations, under a redefini- 
tion of variables, is known to describe two decoupled Riemann equations. This 
case is slightly tricky, since in this case, both the dynamical variables, u, v have 
the same dimension. We describe the Lax description for this separately, which 
does not fall into the above category. Consider the Lax function 

2 1 

L = + -u + -v^p-'^ (30) 

3 6 

The presence of the p~^ gives it a different character from the earlier construc- 
tion. However, it is straightforward to check that the Lax equation 




leads to 



= n w 



X «+^5n+l(f"7r„^), <5n(w7r„2)}^ 



n+1 



■w. 



- -5„+i(w7r„ )g'^{w Ti^ ) 
n+1 



(28) 
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Uf = —UUx — VVx 

vt = -(uw)x (32) 

which are the polytropic gas equations for 7 = 3. The two sets of conserved 

charges, in this case, are obtained from the two possible ways of defining the 
residues of L"+2 around p = 0, 00 respectively. We also note here that, for the 
case 7 = 2 (n = 1), the polytropic gas equations only have a nonstandard Lax 
description. 



5 Conserved charges: 

A Lax description of an intcgrable system has the advantage that the conserved 
charges can be obtained from residues of fractional powers of the Lax function. 
The polytropic gas, on the other hand, is known to have two infinite sets of 
conserved charges [9]. In the non-standard Lax description of the polytropic 
gas, it is known that the two sets of charges arise naturally from calculating the 
residues around two distinct points. In a standard Lax description, however, 
the residues are unique. It is interesting, therefore, to see how the two sets 
of conserved charges will arise in this description. Let us note that a simple 
dimensional analysis shows that we can assign the dimensions [v] = 2, [u] = n 
so that [Ln] = n. From the explicit forms of the two sets of known conserved 
charges [9], it is clear that if they are obtained from our Lax function at all, 
they should arise from the fractional powers and k — ^ respectively, where 
fc = 0, 1, • • • for the first set while k = 1, 2, • • • for the second set. 

Let Hn^k = Resfifi^") and _F„_fe(u, itj, z) = (g„(?i; z^) + ?t. u- z")*^"*" " . Let 



Pn,k{u,'w) be the coefficient of 2:*^"+^ in the Taylor expansion of Fn,k around 
2; = 0. Then: 

ifn,fc = n-(*'+*)p„,fc(u,w;) (33) 



Since 



it follows that 



(w, w,z) = {kn + 1) z" F„,fc_i (u, w, z) (34) 



dpn,k( A / if fc = 



du \ {kn + l)pn,k-i{u,w) ifA;>0 

and, after repeated integration: 

k ^ / k \ 
Pn,/c(u,'w) = ^ — ( ln+l \ p„^k-m{0,w)u"' (36) 

m=0 ■ \l=k-m+l / 
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Using the fact that Pn,2m(0, w) = r„,2m in conjunction with Proposition 
3 we get: 

L2J fT\mn ( ^ \ 

Pn,k{u,w) = ^ I I n /n + 1 y,rnu+X^k-1m (37) 

m=0 ^ \;=TO+1 / 

and therefore 

s=0 m=o ^ '1=0 

(38) 

Similarly, if Hn,k = Res ^I/n then we can show that 

s=0 m=0 ^ 1=0 

The two sets of conserved charges simply correspond to / Hn,k dx and / Hn,k dx 
and coincide with the known conserved charges constructed earUer from a non- 
standard Lax representation upto an overall normalization. This construction 
shows that the two sets of conserved charges can be obtained from the residues 
of two distinct families of fractional powers of the Lax function in this standard 
description. 



6 Generalization to dispersive cases: 

It is an interesting question to ask which dispersive intcgrablc models reduce 
in the dispersionless limit to the polytropic gas dynamics. It is, of course, 
well known that the two boson equation, in the dispersionless limit, goes to 
the polytropic gas equation for n = 1 (7 = 2) and that, for n = 2 (7 = 3), 
the polytropic gas equation is equivalent to two decoupled Riemann equations 
which can be thought of as the dispersionless limit of the KdV equation. It is, of 
course, clear that there may be several dispersive models whose dispersionless 
limit will give the same equation. However, our interest is to find even one 
family of such models. Surprisingly, beyond n = 4, we have not found any 
dispersive generalization of these systems. 

It is well known that the polytropic gas equations for 7 = 2 (n = 1) can 
be thought of as the dispersionless limit of the two boson equation, which is 
described by the Lax operator 

L = d-u + d-^v (40) 
and the nonstandard Lax equation 

(41) 
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Wc will not discuss this system any further. 

The case n = 2 is tricky even at the dispersionless limit, as we have already 
pointed out. We note that the Riemann equation is well understood to be the 
dispersionless limit of the KdV equation, whose Lax description is one of the 
first to have been studied [14]. However, for n = 2, we have two decoupled 
Riemann equations in the dispersionless limit. The dispersive generalization of 
this system is not at all clear. However, with some work, we have found that 
the Lax operator 

L = d'^ + lu+ ld-\d + ^d-hd-'^v (42) 
2 2 4 

through the nonstandard Lax equation 




leads to 



Ut 



UUr 



VV:r. 



(44) 



With a simple redefinition of variables to u± v, this becomes two decoupled 
KdV equations, which, therefore, provides a trivial dispersive generalization of 
the polytropic gas equations for 7 = 3 (or n = 2). It is worth pointing out here 
that this Lax operator also gives through 



dL 
'di 



>i 



(45) 



leads to 



Ut = -Ur,, Vt = -Vx (46) 

which are the elastic medium equations for 7 = 3. Under this description, they 
do not seem to pick up any dispersive terms. 

For 7 = 4 (n = 3), we expect the dispersive generalization to be related to 
the Boussinesq hierarchy, simply from the counting of dimensions. In fact, if we 
choose 



L = d^ ^-vd + u 



the standard Lax equation 



leads to 



5L 
'di 



(47) 
(48) 
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'^t — g (j^Uxxxx ^'^xxxxx ^'^xxx'^ ^^'^xx'^x ~\~ ^{p'^x^x ~l~ ISuU-x 4v Vx) 



Vt = - {2Uxxx - Vxxxx - {v'^)xx + 4.{uv)x) 



(49) 



These clearly provide a dispersive generalization of the polytropie gas equations 
for 7 = 4. These equations are known to admit the W3 algebra as one of the 
Hamiltonian structures. 

Let us further note that the Lax equation 



dL 
'dt 



(50) 



leads to the set of equations 



2vvx) 



Vt = -Vxx + 2Ux 



(51) 



which provides a dispersive generalization of the elastic medium equations for 
7 = 4. In general, we note that the equations 



dL 
'dt 



(52) 



with k = 3m ± 1 defines consistent equations corresponding to the two different 
hierarchies. 

For 7 = 5 (n = 4), consider the Lax operator [15] 



This leads to consistent equations through 



dL 
'dt 



(53) 



(54) 



for k = 4:171 ± 1. Thus, for k = 5, the dynamical equations turn out to be 



1(12 
32 ^ 



^xxxxx ^'^xxxxxxx ^{'^'^'i^xxxx H" '^'^x'^xxx ^^xx^x ~^ 20(^x^)a 

+40uux - ^{v'^Vxx + vvl)x + 5v'^Ux) 



= 1 {40Uxxx - 'i-SVxxxxx - 15{2vvxx + vl)x + 4Q{uv)x - l^V^Vx) (55) 

With a simple change of variables, u = u — ^v'^,v = v, it is easy to check 
that these eqiiations reduce, in the dispersionless limit, to the polytropie gas 
equations with 7 = 5. Therefore, this model provides a dispersive generalization 
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of these equations. It is worth noting here that this system of equations admits 
as a Hamiltonian structure (after suitable redefinition of fields) the nonlinear 
W algebra, W(2,4), which is uniquely characterized by the presence of a spin 2 
Virasoro field as wcU as a spin 4 primary field. 

Similarly, it can be easily checked that, for /c = 3, the Lax equation leads to 

_3 X 3 

— fJ'xxx ^Vxxxxx ^y^'^xxjx ~^ ^UxV 

5 3 

Vt = -^Vxxx + ^Ux- ^VVx (56) 

which, under the same redefinition of variables, goes over to the elastic medium 
equations for 7 = 5 in the dispersionless limit. This system of equations, there- 
fore, gives a dispersive generalization of both these systems. 

For 7 = 6 (n = 5) as well as 7 = 7 (n = 6), we have explicitly verified that 
there is no standard Lax equation that leads to consistent equations with two 
dynamical variables. For higher values of n (and, therefore, 7), it is unlikely 
that a standard Lax description would lead to consistent equations with two 
dynamical fields, since the number of consistency conditions increases rapidly. 
However, we have not checked this explicitly beyond n = 6. The dispersive 
generalization of the polytropic gas equations for higher values of 7, therefore, 
remains an open question. It is possible that they arise only as nonstandard 
equations or that one may have to introduce additional dynamical variables, 
which, somehow, disappear in the dispersionless limit. 



7 Conclusion: 

We have derived a standard Lax description for the polytropic gas dynamics. 
The Lax function, in this case, is intimately connected with Lucas polynomi- 
als, which are also related to the Fibonacci polynomials. The two infinite sets 
of conserved charges have been obtained from the residues of two distinct se- 
quences of fractional powers of the Lax function. We have shown that the same 
Lax function also provides a standard Lax description for the clastic medium 
equations. In addition, we have presented some results on possible dispersive 
generalizations of such systems. 

This work was supported in part by US Deparment of Energy grant number 
DE-FG-02-91ER40685 and by CNPq-Brasil. 



A Standard Lax description for elastic medium 
equations: 

In this appendix, we show how the same Lax function of (20) leads to a standard 
Lax description for the elastic medium equations [7, 9] . 



(-"-■), 



-1 ^n-l 

7r„ 



n-1 
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^i-l n-l 



ri' 



{gn{w z'^)y 



n-l 



{9n{z)f " (z = u;7r„2) 




= n^-i(n-l)t«,7r2"-4 



(-1) 



n-l 



(n - 1) w''' 



Using these results, the Lax equation 

d Ln i_2_ 1 

— -n - 



dt n-l 
can be written 

dLn dL 



Wt^+Ut = u,—^ + (-1)" w"-' 
ow ow 
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which holds if and only if 

{ Vt = Ux 

\ ut = -t;"-^ Vx 

These are none other than the elastic medium equations. Thus, we have shown 
that the same Lax function also provides a standard Lax description for the 
elastic medium equations. 
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